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We present a detailed study of the BCS pairing transition in a trapped polarized dipelar Fermi 
. . . gas. In the case of a shaUow nearly spherical trap, we find the decrease of the transition temperature 

' as a function of the trap aspect ratio and predict the existence of the optimal trap geometry. The 

, latter corresponds to the highest criticai temperature of the BCS transition for a given number 

. of particles. We also derive the phase diagram for an ultracold trapped dipolar Fermi gases in 

(N . the situation, where the trap frequencies can be of the order of the criticai temperature of the 

^ ^ BCS transition in the homogeneous case, and find the criticai vaine of the dipole-dipole interaction 

. energy, below which the BCS transition ceases to exist. The criticai dipole strength is obtained as a 

' function of the trap aspect ratio. Alternatively, for a given dipole strength there is a criticai vaine 

of the trap anisotropy for the BCS state to appear. The order parameter calculated at criticality, 
exhibits nover non-monotonie behavior resulted from the combined effect of the confìning potential 
and anisotropie character of the interparticle dipole-dipole interation. 
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I. INTRODUCTION 



' One of the most challenging goals of modem atomic, molecular and optical physics is to observe the superfluid 
(BCS) transition [1] in trapped Fermi gases. The possibility of such transition for gases with attractive short range 
Q ' interactions has been predicted in Refs. [2] , and has been a subject of very intensive experimental investigations since 
O then (for the latest experimental results see [3]). In typical experiments evaporative cooling is used to cool fermions. 

However, since the Pauh principle forbids the s-wave scattering for fermions in the same internai state, Fermi-Fermi 
y—{ [ [4] or Fermi-Bose [5] mixtures bave to be used to assure collisional thermahzation of the gas. Such combination of 
^ ■ evaporation and sympathetic cooling allows to reach temperatures T ~ O.lTp, where Tp is the Fermi temperature 
' at which the gas exhibits quantum degeneracy. Unfortunately, criticai temperatures for the BCS transition. Te, are 
predicted to be much smaller than Tp. Nowadays, the standard way employs a Feshbach resonance in order to increase 
the atomic scattering length to larger negative values. Such "resonance superfluidity" should lead to superfluid 
transition at Te ~ O.lTp [6]. In the most promising scenario, one starts with a molecular condensate formed for 
fls > and modifìes towards the negative values [3, 9]. Another way to achieve the BCS regime is to use the 
\ cooling scheme that can overcome the Pauli blocking, such as appropriately designed laser cooling [7]. Yet another 
. promising route is to place the Fermi gas in an optical lattice and enter the "high Te' regime [8]. 
d ' The temperature of the BCS transition in a two-component Fermi gas depends dramatically on the difference of 
^ , concentrations of the two components, which presents another experimental obstacle [9]. This problem, however, is 
I ■ not relevant for a polarized Fermi gas with long-range interactions, such as dipole-dipole ones. The possibility of the 
Cooper pairing has been predicted in Refs. [10] and the criticai temperature (including many-body corrections), as 
well as the order parameter, bave been obtained for a homogeneous gas in Ref. [11]. Possible realizations of dipolar 
Q ' gases include ultracold heteronuclear molecular gases [12], atomic gases in a strong DC electric field [13], atomic gases 
with laser-induced dipolcs [14], or with magnetic dipoles [15]. For dipolar moments d of the order of one Debye and 
densities n of lO^^cm^"^, Te should be in the range of lOOnK, i.e. experimentally feasible. 

Dipole-dipole interaction is not only of long-range, but also anisotropie, i.e. partially attractive and partially 
$H repulsive. Thus, the nature of the interaction for trapped gases may be controlied by the geometry of the trap. For 
a dipolar Bose gas in a cylindrical trap with the axial (radiai) frequency ujz{ujp), there exist a criticai aspect ratio 
A = {ujz/ujpy^'^ , above which the Bose-condcnsed gas collapses if the atom number is too large [16], and below which 
the condensate exhibits the roton-maxon instability [17]. The trap geometry is also expected to control the physics 
of trapped dipolar Fermi gases. So far, however, only partial results were known [18]: analytic corrections to Te in 
"loose" traps, and solution of the case of an infinite "slab" with ojp = 0, and ojz finite. In the latter case there exists 
a criticai frequency above which the superfluid phase does not exist. Very recently, we have reported results for the 
case of a general trap [19], and predicted the criticai dipole strength, below which the BCS pairing transition ceases 
to exist. 

In this paper we present the detailed derivation of the results of Refs. [18] and [19], and study the effect of the trap 
geometry on the BCS transition in trapped dipolar Fermi gases. We first consider the case of a shallow nearly spherical 
trap with the trap frequencis Uz, ujp much smaller than the criticai temperature Te of a spatially homogeneous gas 
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with the density cqual to the maximal density of a trapped gas samplc, LOz,LiJp ^ Te- In this case, the presence of a 
confining potential results in a decrease of the criticai temperature as compared to the spatially homogeneous gas. In 
the case of a strong confining potential, where w^, ujp ^ Te, we calculate the phase diagram in the piane F — A""'^, where 
r = 36n(f/n^ is the strength of the dipole-dipole interaction relative to the chemical potential /x. Below the criticai 
vaine of the interaction, F < Fc, the BCS transition does not take place. Similarly, for a given dipole interaction 
strength there is a criticai vaine of A"""^, above which the BCS state appears. We determine the dependence of Fc 
on A~^, and calculate the order parameter at the criticality. The order parameter exhibits a novel non- monotonie 
behavior in strongly elongated cylindrical traps. 



II. BCS PAIRING IN A DIPOLAR FERMI GAS 

We consider a dipolar Fermi gas polarized along the z-direction in a cylindrical trap. The corresponding Hamiltonian 
reads 



-^ + C/trap(r)-M 



m (1) 



+1 ^ ^ V^t(r)v;t(r')ydip(r - rOV^(r')V^(r) 



where 'tp''{r) and 'tp{r) are the canonical fermionic creation and annihilitation operators of particles with the mass m 
and the dipolar moment d, f/trap(r) = m[u)p{x'^ + y^) + t^lz^] the trapping potential, fi the chemical potential, and 

Vdip(r) = {(Plr^){l - 3z^/r^). 

The BCS pairing is characterized by the order parameter A(ri,r2) = Ydi^i — Y2) (^^(ri)^(r2)^, which attains 

nonzero values below the criticai temperature T^. Just below Tg, the order parameter is a nontrivial solution of the 
BCS gap equation (see. e.g. [1]) 

A(ri,r2) = -ydip(ri -r2) / i4:(ri, r2; rs, r4)A(r3, ri) (2) 

>'r3,r4 

and is an extremum of the functional 

^[^h/ tP?^-^+[ A*(ri,r2)K(ri,r2;r3,r4)A(r3,r4). (3) 

7ri,r2 ^dip(,ri - r2j Jri...r4 

The kernel K in the above expression is 

ii:(ri, r2; r3, r4) = T^C;^ Jri,r3)C?_a,„(r2,r4), (4) 

where 



G^Jr,r') = $^ 



iujn - {e^ - fx) 



is the Matsubara Green function of the Fermi gas in the normal phase. In this formula, tOn = 7rT(2n + 1), n = 
0, ±1, ±2, . . ., ipui^) and are the eigenfunctions and the corresponding eigenvalues of the free particle Hamiltonian 
in the trap, [-fi,2v2/2m + Vtrap(r)]V^^(i") = e^tpt,{r). 

The solution of the general problem given by Eqs. (2)-(4), is very difficult even for numerical methods due to 
the large number of variables, relatively low symmetry of the system, and a long-range character of the interparticle 
interaction. Remarkably, in the case of a shallow nearly spherical trap, the solution can be found analytically, whereas 
the general case can be treated by using a variational approach. 



III. SHALLOW NEARLY SPHERICAL TRAP 



We begin with the case of a weakly deformed spherical trap with the frequencies, which are much less than criticai 
temperature, LVp,LVz <C T^. In the new variables 



R=(ri+r3)/2, r = ri - r3. 
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R' = (ra + Ti)/2, r' = Y2- v^, 

the second terni in Eq. (3) reads 



'R,R',r,r 

where 



/ A*(R + r/2,R' + r72)i^(R,R',r,r')A(R-r/2,R'-r72), (5) 

ÌR,R',r,r' 



(R, R', r, r') = T ^ (R, r)C;_,„ (R', r'). (6) 

The kernel K depcnds on variables R and R' only due to the presence of the trapping potential, but, as a function 
of r and r', the kernel deeays rapidly for r, r' > where = Pf /fnTc — vp /Te with pp — mvp — y/2mfl being the 
Fermi momentum, determines the eharacteristic scale for pairing correlations. Under the condition ijjp,ujz <C Te, the 
gap A(ri, ra) = A(Ri2, ria) is a slowly varying funetion (on the scale ^o) of R12 = (ri + r2)/2 (see the end of this 
Section). At the same time, the Fourier transform of A with respect to ri2 = ri — r2, 



A(Ri2,p)= / A(Ri2,ri2)exp(-ipri2), 

ri2 



(7) 



varies on a scale of the order of the Fermi momentum, p ~ pp ~ hv}/^ , see Ref. [11]. It is therefore convenient to 
write Eq. (5) in the form 



A*(R, + P + |) exp(zqr,-*P(r - r'))X(Rc + ^, Re - y , r, r')A(Re - P - |), (8) 



/ / / 

where Re = (R + R')/2 = (R12 + R34)/2 and Te = R — R', and, keeping in mind that the pairing takes place in the 
centrai region of the gas sample, where ?7trap(R) ^ fJ-, we can expand the order parameter in powers of (r + r')/4. 

The leading term of this expansion is 

/ / / À*(Re,P + §)exp(iqre-iP(r-r'))TVG^„(Re + §,r)G_^„(Re-^,r')À(Re,P-^) 

JRc ^re,r,r' JP,q ^ ^ J, A A 

with 

G„„ (R, P) = ^ (R, r) exp(-zPr). 
In the case cOp,u}z <C Te, the Green function Gt^,^(R, P) can be approximated as 

G„„ (R, P) « -—^—^^--^ = ■^^_(^p2/2m-n(R)) 
with /i(R) = — f/trap(R)- With the help of the formula 

tanh(ep/2Te) 



^ ÌUJ„ -^P + Utrap(R) -ii^n - + C4rap(R') ~ L 2/i ^^j 



2Cp 



tanh(^p/2Te) 1 dt^rM^p/^ .... 

--(t^trap(Rej + t>trap(re/^jj;^- , (llj 

/i Ót,p Z^p 



the integration over Fc in Eq. (9) gives 

^ exp(iqrep)T ^ (Re + re/2, P)C?-a,„ (Re - re/2, -P) 
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, 1,, ^ d \ tanh(|p/27y 1 1 9 tanh($p/2re) 
1 - -[/t.ap(Rc)^ j 5{c0 + -U,.^^i-V^)S{q) — — . 



(12) 



The term containing J7trap(Vq/2) can be neglected because Vq ~ 1/pf, and, therefore, Eq. (9) can finally be written 
in the forni 



(13) 



This exprcssion corresponds to the locai density approximation with n /i(Rc) = M — C^trap(Rc) expanded in powers 
of C^trap(Rc)/M up to the first order. 

The next terni of the expansion of Eq. (8) in powers of (r + r')/4 is quadratic (the linear in (r + r')/4 contribution 
vanishes due to the symmetry of the problem) and has the forni 



Re Jrc JP,q 



A*(Re,P + q/2) 



VjVj + V^VJ - 2V^VJ 



A(Re,P-q/2) exp(iqrc) 



r + r \ / r + r 



exp (-iP(r - r'))K(Rc + rc/2, Re - rc/2, r, r') 



(14) 



where Vi and Vi are the i-th componcnt of the gradient Vr^. acting on the left (on A*) and on the right (on A), 
respectively. After neglecting the (^-dependence of A, the integrations over Te and q are straightforward and Eq. (14) 
can written as 



ViVj- + ViVj- - 2ViVj 



A(Rc,P)) 



Jr,r' 2 



r + r 
2 V 4 



r + r' 



exp (-iP(r - r'))ii:(Rc, Re, r, r'). 



(15) 



One can show with the help of the explicit forni of the Green functions, Eq. (10), that the main contribution from 
the integrals over r and r' is 



1 /x(Re) 



1 



-riin 



tanh(^p(Re)/2Te) 



32 ' mT^ cosh2(^p(Re)/2re) 2^p(Re) 

with n bcing the unit vcctor in the direction of P, n = P/P, and ^p(Rc) = P^/2m — /i(Rc). This expression decays 
exponentially for |5p(Rc)| > Te, and, therefore, can be approximated by the simpler one in integrals over P with a 
slow varying functions of P 

1 /x(Re)7C(3) 



where ({z) is the Riemann zeta-function. With this simpliflcation, Eq. (15) becomes 



^y^^.p(M(Re))^y - 



(A*(Re,P) 



ViVj + ViVj - 2ViVj 



^(Rc,P)) 



P=npjr(Rc) 



(16) 



where up{fi{'Rc)) = mpF^Rc)/^''^^ is the density of states on the locai Fermi surface. 

After combining together Eqs. (13) and (16) and performing the variation with respect to A*(Rc,P), we obtain 
the gap equation in the forni 



A(Re,P) = - y Fdip(P-P') 

7C(3) 



,,tanh(eP'(Re)/2re) 



2^P'(Rc) 



A(Re,P') 



16,2 -^(MRc)) / ^(n'VRj^A(Re,P') 



P'=n'pjr(Re) 



(17) 



The first Une of Eq. (17) coincides with the gap equation in the spatially homogeneous case, see Rcf.[ll], with R,, 
being a parameter. Therefore, foUowing the method developed in Ref.[ll], we find that the order parameter on the 
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locai Fermi surfacc has the forni A(Rc,P = npp{Ilc)) = A{Ilc)(po{n), where <^o(n) = ■\/2sin((7r/2) cos('i?n)) (see Ref. 
[11] for more details) and the function A(Rc) obeys the equation 

A(Re) = -Ao(Rc)ln^^^gpA(Re) - ^Ao(O) J ^^l{n'){n'Vnj'A{K,), (18) 

where C = 0.5772 is the Euler Constant, Ao(R) = (mpi? (R)/27r^)rd, with being the dipole-dipole scattering 
amplitude (F^ = — 8d^/7r in the Born approximation) . With the account of the expressions for the criticai temperature 

ri"^ and explicite form of i^o(n) in the spatially homogeneous case, see Ref. [11], the equation (18) for A(Rc) takes 
the final form 



where /i = /2 = 1 — S/vr^, = 1 + 6/77^. Note, that in obtaining this equation from Eq. (19), we also expand Ao(Rc) 
up to the first ordcr in J7trap(R-c)/A*, similarly te Eq. (13). 

Eq. (19) is equivalent to the Schròdinger equation for a 3D anisotropie harmonic oscillator. As a result, the shift 
in the criticai temperature due to the presence of the trapping potential is given by the lowest eigenvalue and equals 



tP -T, _^^T^°^ _ 1 /7C(3) 



TÌ°) V 487r2 V 2|Ao| 



(20) 



In the considered case iVi/Tc<^ 1, the criticai temperature in the trap is only slightly lower than that in the spatially 
homogeneous case. 

Just below Te, the order parameter has the gaussian form 

A(Re)cxexp(- ^ Rl/2lli), 

i=l,2,3 



where /aì = [vp /ujì)\J ujì/T^"^^ [7C(3)/i/487r2(l + 1/2 |Ao|)] is the characteristic size in the i-th direction. One can 

see, that Zaì ^ R'tfì where Rj)p = vp/^i is the Thomas- Fermi radius of the gas sample in the i-th direction. This 
justifics our assumption that pairing takcs place only in the centrai part of the sample. On the other hand, we have 
Iaì ^ ^0 and, therefore, the gradient expansion of the order parameter in powers of r + rM s legitimate. 

For a given number of particle N = ii^/SuzUJp in the trap with aspect ratio A = y^o^^T^' have Uz = wA^/^ and 
ijp = wA"^/^, where u = [cOzi^py^^, and Eq. (20) becomes 



where F{X) = V7C(3)/487r2[2^1 - S/n'^X-y^ + ^/T+6/^^^X'^/% The plot of the function F{X) is presented on Fig. 
1. 

We see that there exists the optimal value of the trap aspect ration A* = 0.81, which corresponds to the highest 
transition temperature in the trap. The existence of the optimal value for the trap aspect ratio is a result of the 
competition bctwccn the anisotropie character of the dipole-dipole interparticles interaction and the confinement of the 
gas sample in the radiai direction. The formcr becomes predominantly attractive for smaller values of A (cylindrical 
traps) and, therefore, favours the BCS pairing. The latter, on the contrary, acts on the pairing destructively due to 
the size effect and, therefore, less pronounced for larger values of A. 



IV. CRITICAL ASPECT RATIO 



We have seen in the previous section, that the trap geometry has a more pronounced influence on the BCS pairing 
in a dipolar Fermi gas, as compared to a two-component Fermi gas with a short-range interaction. This is due to the 

fact that the states, which form Cooper pairs in a trappcd dipolar Fermi gas, have diffcrcnt quantum numbers n^. 
Therefore, their energies are different, at least by the amount of ujz- When this difference becomes of the order of T^, 



0.5 1 1.5 2 
FIG. 1: The function -F(A) versus the trap aspect ratio A. 



the pairing is obviously impossible. As a result, the superfluid transition in a trapped dipelar Fermi gas can take place 
only in traps with Wz < w^c, where the criticai frequency uizc is found to be Wzc = l-8Tc [18]. For uiz < Wzc, as can be 
sccn from Eq.(21), the confinement in the radiai direction dccrcascs the criticai temperature as well. Therefore, one 
would expect the existence of the criticai aspect ratio Ac such that the pairing is possible only in traps with A < Ac 

In this Section we study the BCS pairing in the case of a cylindrical trap where the trap frcquencics can be of 
the order of the criticai temperature in the spatially homogeneous gas, u!p,LOz ~ Tc°\ but stili much less than the 
chemical potential, LOp,L0z <C /U- In this case, the BCS gap equation (2) can hardly be tractable even numerical without 
additional simplifications. 

As it was shown in Ref. [11], the BCS pairing is dominated by the p-wave scattering with zero projection of the 
angular momentum on the z-axis, = 0, in which the interaction is attractivc. Contributions of higher angular 
momenta, although present due to the long-range character of the dipole-dipole interaction, are numerically small 
(see also Ref. [10]). In the p-wave channels with = ±1 the interaction is repulsive and, therefore, leads only to 
rcnormalizations of a Fermi-liquid type, and will be neglected hcrc. Therefore, for the considered pairing problem we 
can model the dipole-dipole interaction by the foUowing (off-shell) scattering amplitude 

T,{p,p',E)=pzp'zJi{E), (22) 

where p is the incoming momentum, p' the outgoing one, and ji{E) some function of the cnergy E. The amplitude 
Fi describes anisotropie scattering only in the p-wave channel with the projection of the angular momentum Iz = 0. 
The function 71 (E) obeys the equation 

7,(E) - -MB') = /" ^ME) {^^tStTo - T^Wa) ''<^''' 

that foUows from the Lipmann-Schwinger equation for the off-shell scattering amplitude [20]; 71 (£^) is assumed to be 

negative in order to guarantec the BCS pairing. The cut-off paramcter A cnsures the convcrgcncc of the integrai and, 
in fact, can be expressed in terms of the observable scattering data corresponding to on-shell scattering amplitude 
with p = p' and E = p'^/m. It foUows from Eq. (23) that 71 (S) is inversely proportional to E, 71 (i?) = 71 {2mE)~^, 
with some coefficient 71. Therefore, the on-shell amplitude is energy independent, as it should be for low-energy 
scattering on the dipole-dipole potential (see Ref. [21]). 

The coefficient 71 determines the value of the criticai temperature Te of the BCS transition in a spatially homoge- 
neous gas and can be expressed through the dipole moment d using the results of Ref. [11]. In a homogeneous gas, 
the order parameter has the form A(p) = pzAo with some Constant Aq, and the linearized gap equation implies 



71 (m) 



f dp pI 


tanh - 1 


1 (27r)3 2^p 


2Tc 



(24) 



where = p^/2m— /i, and the bare interaction is renormalized in terms of the scattering amplitude with 71 (/j) = ^i/pf 



at the Fermi energy ep = IJi = p%/2m along the hnes of Ref. 
over p, we obtain the equation on Te. 



[22] (pf is the Fermi momentum). After integrating 



3 l7il^'F 



, 2/Z 8 , TT ^ 

In— ln- + C 

Te 3 4 



(25) 
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where i^f = mpp/2TT'^ is the dcnsity of statcs at the Fermi energy. After comparing the solution of Eq. (25) with the 
result of Ref. [11] for T^, we find 71 = -24d77r. 
In the ordinary space, the scattering amphtude Fiis 

ri(r,r',ii;) = dJ{r)d,,6{r')j,{E), (26) 

where r and r' are the relative distances between the two incoming and the two outgoing particles, respectively. 
Therefore, the order parameter in the trapped gas, A(ri,r2) ~ {ìp{ri)ìp{r2)) , has the form 

A(ri,r2)=a,5(r)Ao(R), (27) 
where r = ri — r2 and R = (ri + r2)/2, and the corresponding equation for the criticai temperature is 



Ao(R) 



EM„.„.(R)M„.„.(RO^(^^/^^^ + *^^^(^^/^^^ 



71 (Ai) 

-fj^ f .^"ì.n. ^l/. exp(iq(R-RO) 

J {2nY J (2^)3 2ep + 9V4m ^ ^ ". 



2(6 + 6) 

■ Ao(R'). (28) 



Here 6 = C(^iO) ^ = {nx,nz,nz) are the harmonic oscillator quantum numbers, £^{n) = 
h [cOz{nz + 1/2) + ujp{nx + % + 1)] — and the function Mnin2 (R) is defìned as 

M„,„, (R) = M(^) (^)M(^) (x)M(^) {y) (29) 

with 

MÌ%,{Z) = i [ipnAz)dzipnAz) - ipnMdz^nM] , M^^^M = ^nM^uAv), (30) 

^^(z) being the harmonic oscillator wave functions. 
The gap equation (28) is stili hardly tractable numerically and, hence, we employ additional approximations. We 

assume a largo number of particles such that the chemical potential /i is much larger than the trap frcquencics, 
jjL ^ ijJz,0Jp- Therefore, while calculating the functions Mnin2(R)) we can use the WKB approximation for the wave 
functions ipn of the most important for the BCS pairing states with energies near the Fermi energy ep = H- Another 
simplifìcation is due to the fact that the BCS order parameter Ao(R) varies slowly on an interparticle distance scale 
n~^^^ ~ h/pF, where pF = \/2m^ is now the Fermi momentum in the center of the trap in the Thomas- Fermi 
approximation. As a result, the pairing correlations are pronounced only bctwccii states that are dose in energy. This 
allows to neglect /Am in the denominator of the second term in Eq.(28) together with rapidly oscillating terms in 
the functions Mnin2(R)- We then obtain (see Appendix) 

Mi%,{z) ^ M^^iiz) « (-1)"^ Vi - iz/lzNr\Jr.-i{zM (31) 

and 

MÌ'XM ^ M^nÌÌ^) « (-l)"-i / T^xM, (32) 



where n = \ni — <C N = {ni + n2)/2, = y/2Nh/mji — IqìV^N, whereas U„(2;) = sin((n + 
l)arccos2;)/ sin(arccos z) and T„(a;) = cos(n arccoscc) are the Chebyshev polynomials. The functions M^j^jC-^) and 
Afilins {^) fuUfil the foUowing completeness relations 

E = - (33) 

n>l 

E^nM(f)(x)M(f)(x') = ^ / J {x-x% (34) 

n>0 y^pN ^ 
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with So = 1 and (5„>o = 2, which foUow froni the coniplcteness of the Chebyshev polynomials. 
It is convenient to rewrite Eq. (28) in the foUowing way 



Ao(R) 

71 (m) 



JéR' [e M,,.(R, | ».-.«./2T) + ..nhfc/2r) _ tanh((C.+&)4T) | 



-E/W(R){i=H!|l±|p-^|M„.„,a., 

J dR' [Ki(R, R') + K2{R, R') + iiTaCR, R')] Ao(R'), 



)3 2(p2/2m - A*) 



cxp(iq(R - R')) 



Ao(R') 



(35) 



because the sum + ^2 does not depend on m — n2 and, therefore, with the help of formulae (33) and (34), the 
summation over n can casily bc pcrformcd in kcrncls K2 and K3. On the other hand, the kernel Ki is determined 
entirely by the states near the chemical potential /x. 

The calculation of the sums in the kernel ii'3(R, R') gives 

^ M„,„,(R)-^M„,„,(R') = ^ ^ Mn„(R)-^Mn„(RO 



- ^(-R Ti ) 1 (»^^)^ 

- 6{R-R) ^ 2^(N) ^3 

dpi pI 



/2 - z2 



Jì'\ \^ 1 {nUJz? IzN 



.(R-R')/nt 



2(pV2m-/i(R))' 



where ^(R) = — X^i ^^ÌRÌf^ì ^ind we havc rcplaccd the discrete sums over Ni with the integrals over contineous 
variables Pi = \/2muJiNi. As a result, the kernel K^i^, R') is 



2(j?2/2m - /i(R)) 2(p2/2m - /i) 



5(R - R'). 



(36) 



FoUowing Eq.(23), we scc that the kernel results in the replacement ^ — > /u(R) in the scattering amplitude 71 (/u) 
on the left-hand-side of the gap equation (35). 

We are interested in the criticai value Ac of the aspect ratio, below which the BCS pairing does not take place for 
a given strength of the dipole interaction. Therefore, this vahie corresponds to vanishing criticai temperature. We 
therefore calculate the kernels Ki and K2 in the limit T <C ujì- The summation over ^(N) = (^1) +^i)/2 in the kernel 
K2 is then within the limits — /x < ^(N) < —Wz/'2,, where —oJz/'^ is the upper limit due to the fact that the function 
is nonzero only when \niz — n2z\ > 1- These sums can again be replaced by integrals with the foUowing result 



(^) 



i^2(R,R') « ^pUR)M^) {lii^^ - |(4 - 21n2)| 5(R- R'), 



(37) 



where pf(R) = -\/2m/x(R) is the locai Fermi momentum in the Thomas- Fermi approximation. 

In order to calculate the kernel Ki in the limt T loì we note, that nonzero contributions to the sums over ni and 
n2 originate only from the region W2/2 < |^(N)| < e{a)/2, where tiì = \nu — n2i\ and £(n) = fi^^Wini. As a result, 
the kernel Ki (R, R') can be written as 



^^i(R,R') 



E 



n^>0;n^,ny>0 tj^/2<?(N)<e(n)/2 ' 



Mn„(R)Mn„(R') 



E ^njuj^' ^^<5(C-^(N))Mn„(R)Mn„(R'). 
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For a gcncric cylindrical trap, /(^, n) = X]n^(C ~ C(N))MNn(R)A^Nn(R-') is a smooth function of ( (this is 
equivalent to neglecting the so-callcd shell efFect). At the same timo, only small n with £(n) <C /x contributes to the 
kernel Ki. We therefore can replace /(C,n) with /(0,n). After replacing the sums over Ni, i = x,y,z, with the 
integrals over a = ujpNx/ (3 = ijOpNy/ fi, and C = '^z^zjii-, we obtain 

i^i(R,R') « - '5„„'^„„ln^l]<5(^(N))MNn(R)MNn(R') 

where Nx = afi/uip, Ny = (ìji/uip, and Nz = QfijLOz- 

We write the order parameter in the forni Ao(R) = Ao(zR^p, xFÓfp, ylìi^p) = Ao(r), where R^^p = pp/rniJi is 
the Thomas- Fermi radius of the gas cloud in the i-direction and |a;| , |?y| , |z| < 1 are dimcnsionless variables. After 
combinning together Eqs.(35),(38),(37), and (36), the gap equation in limit T <C reads 

^(1 - r')Ao{r) = (1 - r2)3/2 L ^ - |(4 - In 4)1 Ao(r) - ^ /' K{r, r')Ao(r'), (39) 
where T = up, /x(r) = M - 14rap(r), and 



Snjn, -4- I dadpdCS{l -a- 13- C)MN„(R)MNn(R'), (38) 



-^(r, r') = V 5njny + — ("a; + %) 



X /" / da I dl3 S{l-C-a- P) 
Jm, Jm^ Jm„ 



4 V(C-.^)(C-.-^) .. . ^ . - . 
V c u„._i(-^)u„._i(^) 



^2 / / r,N/ , tT "J=(~7^) "»(~7^) 

TT"^ y'(a - x2)(a - a;'2) V" V" 
4 1 y y' 

with Ms = max(s^, s'^) for s = x,y, z. 

Ecfore solving the above equation numerically, let us analyse the behaviour of its solutions near the edge of the 
gas sample, r — > 1. In this region max(x^,a;'^) = x^, in&x{y^,y''^) = y^, max{z^ , z'^) — z^, and it is convenient to 
introduce the new variables a = a — x^, b = f3 — y"^, and c = ( ~ z"^. The delta function in Eq. (40) then reads 
(5(1 — — a — b— e) and, therefore, only small a,b,c < 1 — ^ contribute to the integrai. This allows as to write 
the integrai in the last term in Eq. (39) in the form 



i^(r,r')Ao(r') « f^J ^ J2 ^nM^H^z + ^irix+riy)] dadbdcJ -^S{1 - ~ a ~ b ~ c) 



nz>0;nx',ny>0 



* I ^'■y(Y^^72y(Yz^u„,-i(i)Un.-i(^0T^^ 

(1_ ^2)3/2512 5„^5„^ln[n, + ^^(n,+n^)]U„_i(l)T„,(l)T„^(l) 

OTT (jJz 
nz>0;nx',ny>0 



* / V (1 - y,,^ Vn^-i{z')Tr,Ax')Tr,Ay')M^^',yy',zz'), (4i) 

and simplc analysis of Eq. (39) gives the asymptotics Ao(r) ^ (1 — r^y/^ for r — > 1. Eq. (39) in the limit r — > and 
Eq. (41) can be thus used as a consistency test for numerical solutions of Eq. (39). 
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Our numcrical approach to Eq.(39) is based on the observation that it can be considered as the equation on extrema 
of the quadratic form 



F[Ao] = ^ Ao{r)[L{r)S{r - r') - i^(r,r')]Ao(r') 



with L(r) = 3(1 - r2)/r - (1 - r2)3/2{ln([2/z(r)/w^] - 2(4 - ln4)/3}. Taking into account the asymptotics of the 
Solutions of Eq.(39), we will use the ansatz 

Ao(r) = (1 - r2)V2 ^ Cm.m,VmAz^)Tm,{x^ + y^) (42) 
and write the functional F[Ao] as a quadratic form of unknown coefficients Cm^^mp 

-^[^] — ^ ^ -^rriz ,'>Tip,nz ,npC-mz'>Tip^nznp 

nz^np,mz^mp>0 

with 



Mm.,m,,n.,n, = ^ / dr | ^ii^ - (1 - r^)^/^ [in ^ (1 - .2) - ^ (4 - 2 In 2)] | 

^njny H^z + —{nx+ %)] 



/' 

Jo 



dadpdC S{l-a-p-OCVr" '{a,p,C)VA" '{a,P,0, (43) 

where the functions V are defìned as: 



J-l TT ^/l-X^ J-l TT _ y2 TT 



z^Vn^-i{z)UnACz')TnAax' + f3y') (44) 



An extremum of this form at Te corresponds to the eigenvector of the matrix Mm^,mp,n^,np with a zero eigenvalue. 

The condition that the matrix M has a zero eigenvahie imposes a constraint on the interaction paramcter T and 
the trap frequencies Ui and, therefore, allows to determino the criticai aspcct ratio. Indeed, the paramcter T can be 
written as F = 36n{0)d^ /nfi, where n(0) = {2mii)^^'^ /dn^fi^ is the gas dcnsity in the center of the trap, and, hence, 
for a given dipole moment d, T depends only on the chemical potential ^. On the other hand, the chemical potential 
H and the total number of particles in the trap N are also related, 3A^ = /x^/wzo;^. As a result, for a given F and 
A'', the product of the trap frequencies u^uìp is completely determined, and the only free parameter is the trap aspect 
ratio A = (wzjijJpj^l'^ ■ We may thus determine its criticai value Ac from the condition that the lowest eigenvalue is 
zero. 

The calculation of the matrix elements Mm^,mp,,i^,np is naturally divided into two parts (see Eq. (39)). The 
contribution with the locai kernel I/(r) is a three-dimensional integrai that can easily be computed using, for instance, 
the Monte Carlo integration routine, such as the VEGAS algorithm from the GSL library [23]. The part with the 
non-local kernel K[v^ r') is a triple sum, which elements are eight-dimensional integrals. The straightforward approach 
based on the same numerical method fails in this case because it is too time consruning. To overcome this problem, 
we calculate the functions ^^'""''(q:,/^, C) in the following way. We integrate numerically over r' for fìxed a, /3 and 
n's (the vaine of C is thcn fixcd by the J-function) and use these data for a two-dimcnsional splinc intcrpolation of 
the integrand for the last integrations over a and /?. In this way, the time required to compute the matrix elements 
Mjn:,,mp,n:,,np reduces to about seventy two hours on a standard workstation. This procedure gives convergent results 
with rcspcct to the highcst powcrs of the polynomials Mc_ and Mp in our ansatz, Eq. (42). We also chcckcd that our 
Solutions do not depend on the number of shots in the Monte Carlo algorithm and on the mimbcr of points chosen 
for intcrpolation. They are also a proper asymptotic behavior for r ^ 0. 

The results of our calculations are presented in 2 figures. Fig.2 shows the desired relation between the interation 
strength F and the inverse criticai aspect ratio A~^. The three curves correspond to three different numbers of 
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FIG. 2: The criticai lines Fc versus the inverse aspect ratio A ^ for different numbers of particles. The BCS pairing takes place 
above the depicted curves. 

particles. As it could be expected, for larger number of particles, the criticai aspect ratio Ac is smaller because the 
interaction is stronger. For a pancake trap, A~^ < 1, the interaction is predominantly repulsive, and higher values 
of r for a fixcd A are rcquircd to achicvc the BCS transition. On the otlicr band, for a cigar trap, A""'^ > 1 , the 
interaction is predominantly attractive and the BCS transition takes place at smaller values of T. The existence of 
the criticai interaction strength (for a given value of A) is a result of the discreteness of the spectrum in the trap and 
of the spccific structurc of the ordcr paramctcr p^). The lattcr allows pairing only bctwccn particles in the statcs, 
where quantum numbers differ by an odd number (intershell pairing). Therefore, the pairing correlations bave to 
be strong enough to overcome the corresponding energy difference. 

Fig. 3 shows the order parameter Ao(r) for the cigar trap with A"""^ = 2.2. The order parameter is a non-monotonie 
function of the distancc from the trap center, in contrast to the case of the BCS order parameter in a two component 
Fermi gas with short range interactions [24]. This effect persists, although being less pronounced. for the case of a 
pancake geometry. In the axial direction, the order parameter A{z,p = 0) develops a minimum at p < 1, whereas 
in the radiai direction A{z = 0, p) becomes negative in the outer part of the cloud. This completely new behavior, 
originating from the anisotropy of the interpaticle interaction, can bave profound consequences for the form and 
spectrum of the elementary excitations. 



Wc bave prcscnted a dctailed thcory and analyzcd the influence of a trapping potential on the BCS pairing in a 
dipolar single-component Fermi gas. We bave determined the phase diagram for trappcd dipolar Fermi gases in the 
F — A~^ piane, where F measures the dipole strength and A is the trap aspect ratio. The BCS transition at finite 
temperature T is possible iff F > Fc(A). We bave calculated the criticai line Fc(A), and the order parameter at 
criticality. 



We acknowledge support from the DFG, the RTN Cold Quantum Gases, ESF PESC BEC2000+, the Russian 
Foundation for Basic Research, QUDEDIS, INTAS, and from the Alexander von Humboldt Foundation. 



Here we present the derivation of the WKB-expressions for the functions Mnin2{x) and Mnln2{z), Eqs. (32) 
and (31). Our starting point is the WKB-eigenfunctions (Pn{x) of a one- dimensionai harmonic oscillator with the 
Hamiltonian 
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FIG. 3: The order parameter for the aspect ratio A = 0.45 (cigarshape trap). The sohd hne shows Ao(z, p = 0), and the dashed 
Une corresponds to Ao(2 = 0, p). 



with the eigenvalues i?„ — hjj{n+ 1/2), 



' 2muj 



cos[$„(a;) -7r/4], 



(45) 



where 



and 



Pn{x) = \j2m{En - mup-x^ 11) 



*"(a;) = - / Pn{x')dx' = -n + - pn{x')dx' . 

As it was already mentioned, the functions M^ln^ (x) and Mn%.2 (z) with ni and n2 close to each other (and both 
are much larger than unity) are the most important. We therefore introduce N — (ni +n2)/2 and n — ni — where 
n ^ N, and write the expression for Ali^li^^x) as (see Eq. (30)) 

^"^^ {cos[$Ar+„/2(a;) - $Ar_„/2(a;)] + cos[$7v+«/2(a^) + '^N-n/2{x) - 7r/2]} 



ttPnÌx) 2 



■kPn{x) 



{cOs[<I>A,+„/2(x) - <I'Ar-„/2(a;)]} , 



where we neglect the rapidly oscillating contribution cos[$7v+n/2(2;) + ^^-«72(2^) ^ I^i the case n <C iV, the phase 
difference ^N+n/2Ìx) — ^N-n/2(x) can be simpUfied as 

^N+n/2{x) - ^N^n/2{x) = 2"+^/ \PN+n/2{x') - PN-n/2Ìx')]dx' 

^ TT r dmix^ , 



TT — arccos 



(^xy/2Nmuj/h 



(46) 



As a result, we obtain the expression 



Mi^i (x) « (-1)" , cos ["n arccos (xJ2Nmùjlh] 

itpn{x) i \ J 

= (-1)"T„ (.Tv/2A^mcc;/?i) , 
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which corrcsponds to Eq. (32) in the dimensionai nnits. 

The expression (30) for the function Mn^n^iz) involves the first derivative of the wave function. According to 
the general rules of the WKB approximation (see, e.g., [25]), this derivative should act only on the rapidly varying 
trigonometrie factor in Eq. (45). Therefore, Mnin2{z) can be written in the form 

-^niLl^) = ^nÌÌ^) = \ [VN+n/2{z)dz^N-n/2{z) - ^N-n/2{z)d^^PN+n/2{z)\ 
~ -^ {^™[^N+n/2{z)-^N-n/2{z)]] , 

where again we neglect the rapidly oscillating contribution. The application of Eq. (46) gives 

MifnA^) « (-l)"^sin narccos (^x^/2NmoJ/flj 

and, assuming n > 0, we obtain 



Mill,{z) « (-l)"^y 1 - (^x^/2Nmw/h) U„_i (^x^/2Nmu/h) . 

Note also, that one obtains the same result if, instead of differentiating the WKB wave function, Eq. (45), one uses 
the well-known relations between the wave functions of a harmonic oscillator and its derivatives. 
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